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ABS TRACT 


An investigation is made of the microprogramming approach to 
logical design of a digital machine. The digital processor designed 
performs general matrix manipulation and in particular is adapted to 
the requirements for implementing a Kalman-Weiner filter in a Track- 
While-Scan radar system. 

A timing and data-flow analysis is made and micro-programs written 
to implement the required macros for the selected benchmark applica- 
tion. A detailed comparison is made of the operation of the micro- 
programmed processor to that of a current general purpose avionics 
type computer of similar main memory cycle time. Some conclusions 
are made concerning the optimizing of various parameters of a micro- 


programmed computer design. 
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Symbol 


LIST OF SYMBOLS 


Definition 
A Register (16 bits) 


The contents of the lower half, (Ag...A,_,,), of the 
: 9 16° 
A Register 


The "a" field of a main instruction word containing 
a five bit opcode 


The "b'' field of a main instruction word containing 
the first eight bit operand address 


Main Memory Address Register (8 bits) 


The "c'' field of a main instruction word containing 
the second eight bit operand address 


The "d'' field of a main instruction word containing 
the eight bit result address 


Microinstruction Decoding Register (16 bits) 


The "f'"' field of a microinstruction word containing 
the five bit opcode 


Control Memory Address Register (8 bits) 
I Counter Register (16 bits) 
Initialization Register for I (16 bits) 


The "i"' field of a main instruction word containing 
a three bit index 


J Counter Register (16 bits) 
Initialization Register for J (16 bits) 


The "j" field of a microinstruction word containing 
a three bit y field modifier 


K Counter Register (16 bits) 
Initialization Register for K (16 bits) 


The "a'' field of a microinstruction word containing 
eleven bits used with concurrent operation opcodes 


Symbol 


Z 
Zo 


CY) 
M<x) 


(N{Z)) 
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Definition 


The elements of the k field 


Main Memory (256 x 16) 

The column dimension of a matrix 

Control Memory (256 x 16) 

The row dimension of a matrix 

Q Register (16 bits) 

Main Instruction Decoding Register (16 bits) 

S Register (16 bits) 

A cell in Control Memory used for temporary storage 
A cell in Control Memory used for temporary storage 
X Register (8 bits) 

Initialization Register for X (16 bits) 

Y Register (8 bits) 

Initialization Register for Y (16 bits) 


The "y'' field of a microinstruction word containing 
an eight bit operand address 


Z Register (8 bits) 
Initialization Register for Z (16 bits) 
The contents of the Y Register 


The Main Memory cell addressed by the contents of 
the X Register 


The contents of the Control Memory cell addressed by 
the contents of the Z Register 


Indicates a transfer operation 


Symbol Definition 


(Az. +-Ay) The contents of the ith through the jth bit position 
of the A register. Register bit positions are nun- 
bered from left to right starting with one 

Dero ed Operation q is performed if p is true 

58— A The octal number 58 is put in the A Register 


(A) ——9302 The contents of the A register are transferred to 
cell #302 in N 


1. Introduction 

The great majority of general purpose digital computers are designed 
around the concept of sequential performance of stored instructions. 
These instructions may be as simple as merely transferring a quantity 
from one place to another or as complicated as multiplying two quanti- 
ties together. The process of performing one instruction is usually 
subdivided into two parts or cycles. The first cycle is called the 
"Read Next Instruction" (RNI) cycle and refers to the transfer of an 
instruction word from memory to a special instruction register. The 
second cycle is the "Execution" (EXEC) cycle and is concerned with the 
interpretation of the instruction word and execution of the specified 
operation. 

Any computer operation can be reduced to a series of data trans- 
fers between various registers and/or memory locations. These trans- 
fers are controlled by opening and closing gates in the data transfer 
paths. Thus the execution of a specific instruction requires the 
generation of a specific sequence of transfer command signals. This 
process iS represented in Figure 1 where the inputs are the code for 
the desired operation (opcode) and the address(es) of the operand(s ) 
which are all contained in the instruction word. Of particular inter- 
est here is the method of implementing the "transfer command generator". 

The Conventional Design Approach 

The conventional manner of generating these signals is through 
the use of hardwired logic circuitry consisting of AND-OR-INVERT, 
NOR-NAND, or equivalent logic and FLIP-FLOPS. (See Figure 2). The 
instruction word is transferred to the instruction register where 


the opcode and address portions are decoded by separate decoding 


matrices. A clock, which controls the timing of the basic computer, 
generates timing pulses which cause the timing and sequencing control 

to step through a series of states at precise time intervals. Each 
state generates a signal which is combined with the outputs of the de- 
coding matrices in a control matrix. The control matrix in turn 
generates the proper gate control signals to execute each data transfer. 

As the number of memory locations, registers, and opcodes in- 
creases, this circuitry rapidly becomes very complex and usually is a 
major portion of a computer's hardware. Another aspect of the usual 
computer design concept is that the computer designer determines the 
manner in which a given instruction is executed and this, being in- 
variant once the computer is built, places a limit on the flexibility 
of the machine. Usually, to prevent this from being a Serious limita- 
tion, the instruction set is carefully constructed to meet the antici- 
pated needs of the user, but this of course, results in an even more 
complex "generator". One way of overcoming these disadvantages is by 
use of 'Microprogramming". 

The Microprogramming Approach to Computer Design 

The microprogramming concept is not new, but until recently, it 
had not been utilized in commercially available, general-purpose com- 
puters. 

Basically, microprogramming is another method of implementing the 
"transfer command generator". A definition which reflects the most 
characteristic feature of modern microprogramming design practice is 
as follows: 

"Microprogramming extends the concept of a stored-program digital 


machine to include a second, lower or more basic level of stored in- 
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structions for purposes of either flexibility of function, economy of 
computer logical circuitry, or a combination of both". 

This definition implies the use of a second memory containing 
stored-instructions each of which commands a particular data transfer. 
Thus by entering this second memory or control memory, as it is usually 
called, at the proper point and executing a given Sequence of these 
second-level or microinstructions, the necessary series of transfer 
commands can be generated. The main opcodes are constructed to point 
to the proper entry point in control memory. 

The general configuration for microprogramming is shown in Figure 
3. The main opcode goes to the control memory address register and 
initiates the proper entry into the Series of microinstructions stored 
Ln Acont no L memory . The microinstruction is transferred to a micro- 
instruction register. From there, it is applied to a decoding matrix, 
along with the operand address and timing pulses. The output of this 
matrix is the necessary transfer command signals. Each microinstruc- 
tion also specifies the next microinstruction of the desired sequence. 

This, in general, is a brief outline and comparison of the con- 
ventional versus microprogramming approach to computer design. It is 
clear that the microprogramming idea hinges on the availability of a 
Suitable efficient and economical control memory. The increasing 
capability of such small, fast "scratchpad"' memories is evident in a 
number of recent reported designs. [4| The digital machine logic de- 
Signer is presented with the potential for significant improvements in 
machine performance. It is the purpose of this thesis to explore and 
evaluate the microprogramming design approach. To this end the re- 


quirements of a typical real-time digital process are used to establish 


ses 


the parameters of a machine. The design is then carried out, using 
principles of microprogramming. Programs are written for the specified 
application and the overall machine performance compared with that of 

a representative current military real-time control computer of similar 


arithmetic and main memory speed. 
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ae Objective 


In today's world of high speed aircraft and missiles, there is an 
ever increasing need for more complex control systems. These systems 
must be capable of observing and predicting many state variables and 
outputing accurately and quickly the many control signals required. 
The digital computer in conjunction with sampled-data control theory 
is helping to supply this need. Also with the advent of microelectron- 
ics these computers can be made small enough to make them feasible for 
aerospace applications. Thus there is a requirement for a computer to 
do a specialized job. 

A common question that often arises at this point is one of "gen- 
eral purpose" computer versus "special purpose" computer. By "general 
purpose" is meant a computer capable of general scientific computation 
Such aS is usually found in industry, universities, and laboratories. 
"Special purpose" denotes a computer built to handle a specific type 
of problem and therefore is generally very limited as to general- 
purpose capabilities. The intuitive answer would be to choose a 
"special-purpose'' computer for a specialized job. However, this may 
not be the right answer. 

A general-purpose computer can usually be programmed to perform 
the function of any special-purpose computer and therefore is adaptable 
to a variety of special-purpose applications. Also the GP computer is 
more likely to be commercially available, thus saving the user the 
R & D costs associated with producing a new machine of limited use. 

On the other hand, the special-purpose idea is not without its 
advantages. Since it must perform only specified tasks, it can be 


designed to do them with maximum efficiency, and the hardware can be 
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minimized which will give the SP computer a definite advantage with re- 
gard to size, weight, and reliability. It is thought that the micro- 
programming approach to computer design could produce a flexible com- 
puter in a small enough package to be a useful compromise solution and 
it is the objective of this thesis to investigate this hypothesis. 

In order to make a quantitative analysis, it is necessary to as- 
sume a specific situation, assign realistic parameters, and evaluate 
the results. The application of a computer as a digital filter ina 
Track-While-Scan airborne radar system was chosen as the specific 
application. In this case, the range to a target is sampled period- 
ically in the presence of noise. It is desired to filter these range 
measurements or Samples to produce a more accurate estimate of the tar- 
get's actual range at any given sampling instant. R.E. Kalman developed 
a set of recursive equations based on sampled-data theory, commonly 
known as the "Kalman Filter", for this purpose. A block diagram of 
this filter is shown in Figure 4, where the single lines indicate 
scalar quantities, and the double lines indicate vector quantities. 

The actual equations are shown in detail in Appendix VII. A fourth- 
order system was assumed with only one state variable observable. 
Therefore the filter input is a Scalar whereas the output is a vector 
consisting of the predicted values of all four state variables. It 
will be further noted that the filter equations are formulated in 
Matrix notation and it is this fact that led to the design of a com- 


puter which specialized in matrix arithmetic. 
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3. Computer Design 
Hardware 

It is appropriate to begin by considering the basic word structure, 
since many computer parameters are directly related to this character- 
istic. A computer word, composed of a sequence of binary digits (bits), 
generally is interpreted in one of two ways. Either the entire word 
represents data (a number) in binary or it represents a computer in- 
struction. In the latter case the word is generally divided into at 
least two parts, one being the code indicating the operation to be per- 
formed (opcode), and the other giving the address in memory of the 
operand or the operand itself. Thus the first two decisions required 
are 1) how long should the word be, and 2) how to subdivide the word 
when using it aS an instruction. In this case these decisions had to 
be arbitrary, but an attempt was made to make them realistic. 

It was assumed that this computer would be used as a digital fil- 
ter and predictor in a Track-While-Scan radar system and that the 
maximum data quantity would be 1000 miles of range to an accuracy of 
-l miles. To represent this quantity in binary would require 14 bits. 
Allowing one bit for addition overflow and one bit for sign, a total 
word length of 16 bits was chosen. 

The second decision regarding how to subdivide the word is more 
complicated due to its interaction with other facets of the computer 
structure. It is generally desirable to allow any memory location to 
contain an operand. Therefore the number of bits alloted for operand 
address immediately determines the size of the memory or vice versa, 
using the relation that n bits can specify 2" locations. In a similar 


fashion, the number of bits alloted to opcodes immediately prescribes 
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the maximum possible number of operations available. Another side to 
this problem peculiar to the microprogramming concept is the need for 
compatibility with the separate control memory which contains the micro- 
programs. Despite the various dependencies on this decision mentioned 
above, it was actually a software consideration which really influenced 
the final choice. 

It seemed desirable to be able to specify a complete operation in 
one instruction. In other words, an instruction would give the opcode, 
the locations of the two operands (in the more general case), and the 
location of the result. To do this it was decided to subdivide the 
word into 2 equal parts of 8 bits each and actually use 2 successive 
words for each instruction. The upper portion (first half) of the 
first word would contain the opcode, and the lower portion would con- 
tain one operand location. The second word would contain the second 
operand location, if any, and the location of the result in the upper 
and lower sections respectively. An immediate consequence of this de- 
cision is that the main memory is limited to 256 words which may be 
too small, but it was decided to evaluate this at the conclusion of 
the project. 

At this point it is necessary to look at some of the special 
features required to implement the microprogramming concept. The most 
important one is the Control Memory (designated N) which contains the 
microprograms. This memory is separate from the Main Memory (designa- 
ted M) and is usually an order of magnitude faster than M. The reason 
for this latter characteristic will be explained later. The micropro- 
grams appear much like subroutines. They are in effect "called" by 


the opcodes in the main program and are given the operand and result 
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locations aS arguments. The microprograms in turn consist of micro- 
instructions which are executed sequentially. They must be written to 
completely perform the desired operation including the considerable 
bookkeeping inherent in matrix arithmetic. 

Of course the question of how big to make the Control Memory has 
to be decided and here again 256 words were chosen. This means that 
an 8 bit address is required as was true for the Main Memory. This 
latter fact was not a coincidence by any means, but was done deliberately 
for a reason which will be explained. 

The next step in the design process was to look in detail at the 
type of arithmetic operations to be performed. It was decided to con- 
sider the following five matrix manipulations: 

1) Matrix Addition (MATAD) 

2) Matrix Subtraction (MATSB) 

3) Matrix Transposition (MATPS) 
4) Matrix Multiplication (MATML) 
5) Matrix Inversion (MATIN) 

Matrix Multiplication was selected in particular as one of the 
more complicated yet typical matrix operations. It is necesSary to 
store the elements of a matrix in memory in a specified sequence. This 
author chose to store them by rows reading the elements from left to 
right in conventional matrix notation. In the two cells immediately 
preceding the elements, the number of rows, n, and the number of col- 
umns, m, of the matrix would be stored in that order. Thus the matrix 


A would be stored in memory from location L onward, as follows: 
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Memory Location Contents of cell 


L n 
L +1 m 
L + 2 Ay 
L + 3 Aj» 
L +m Ain 
+ + j 
L m Ay 
L + (mxn)+#+il 1a 


Matrix Multiplication (A x B = C), is performed by means of the follow- 


ing algorithm. 


h,=n 
ae A_B ee Ee ee 8) 
= k=1 ik’ kj Dale 25 is Mp 


It was noted that at any given step in the above computation it is 
often the case that different elements of three matrices are refer- 
enced. Therefore to facilitate the acquisition of these elements as 
needed, it was decided to have three address registers, one for each 
matrix. Further, these registers would be able to address any location 
in the Main Memory or the Control Memory; hence the compatibility in 
Size between the two memories. Also the three address registers (desig- 
nated X, Y, and Z) would have the capability of being incremented or 
decremented in one step to facilitate reiterative steps being performed 
on successive matrix elements. 

One final decision required to complete the specification of the 
Control Memory is the word length. Obviously this would have to be 


something greater than eight bits since the Gontrol Memory contains 
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256 words. Although there is no definite reason, sixteen bits was chosen 
to permit temporary storage of little used microprograms in the Main 
Memory if the Control Memory proved to be too small. 

Referring again to the matrix multiplication process, it was found 
that not only was it necessary to reference the elements of a matrix 
sequentially, but often this process had to be repeated more than once. 
To facilitate the reinitialization of an address register, three cells 
in the Control Memory (designated X,, Yo: and Z,) were set aside to 
hold the initial address of each address register and have the built-in 
capability of transferring this quantity to their respective address 
registers in one step. Of course, some type of counter is always essen- 
tial in this kind of operation, so three cells in the Control Memory 
(designated I, J, and K) were set aside as counters with a one-step 
increment or decrement capability. The reinitialization problem occurs 


here also; so In2 J and K, were designated with the one-step transfer 


O 
capability Similar to X,, Yo, and Z,- A sixteen bit arithmetic register 
(A register) was necessary for basic addition and subtraction and to 
provide a built-in multiply and divide capability, two other sixteen 
bit registers (Q and S) were required. As is the common practice, the 
A and Q registers could be connected together to form a double length, 
AQ, register with a long right or long left shift capability for use 
in the multiplication and division process. 

Finally two eight bit registers, C and H, to act as basic memory 
address registers for M and N respectively were provided, along with 
a sixteen bit F register which is used to decode the microinstructions. 


This completed the basic hardware characteristics of the computer and 


the resulting block diagram with data transfer paths is shown in Appen- 
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dix I. The single-ended arrows indicate unidirectional data flow, 
while the double-ended arrows indicate a bidirectional data flow, not 
necessarily over the same wires, however. 

Instruction Format and Operation Codes 

The eight bits available for operation codes in a main program 
instruction word were subdivided into two parts. The first part con- 
sists of five bits and actually addresses the particular cell, in the 
First thirty-two cells of Control Memory, which acts as the entry point 
for the particular microprogram required. The second group of three 
bits acts as an index group. It was found in studying the microprogram 
requirements of MATAD and MATSB that these two operations required 
identical "bookkeeping" and only the actual arithmetic operation was 
different. Therefore the same microprograms could be used provided 
there was a microinstruction within the microprogram to sense the ap- 
propriate arithmetic operation. This microinstruction references the 
three bit index of the main instruction to make this decision. To the 
programmer, MATAD and MATSB would each have a specific eight bit op- 
code, but actually only the last three bits would be different. Thus, 
in general, the use of the index permits minor branching within a given 
microprogram. Only the opcodes for the five basic matrix operations 
were established and they are listed in the Appendix II. Obviously 
there iS a capability for a great many more. 

There were several possibilities considered regarding the format 
of a microinstruction and the sequencing of these instructions. One, 
in particular, was to incorporate into an instruction the address of 
the next instruction to be performed, thus minimizing the number of 


times a particular instruction appeared in the Control Memory. This 
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idea appeared to require too complex a coding scheme to achieve the de- 
sired flexibility, but an evaluation of the idea will be discussed in 
the Conclusions Section. The system actually implemented was the basic 
two cycles used in most computers today. The first cycle is usually 
referred to as the ''Read Next Instruction" (RNI) cycle and does exactly 
what the name implies. It transfers the next instruction from memory 
to a special decoding register. In this particular instance this would 
be a transfer from N to F. The second cycle is called the "Execution" 
(EXEC) cycle and includes the decoding of the instruction word to de- 
termine the required operation and the location of the operand, if 

any, aS well as the performance of the actual operation. A particular 
characteristic of this mode of operation is the sequential execution 
of instructions stored in memory. This may well have been expected 
from. the use of the term “Read Next> Instruction’. “Avfurther implica— 
tion here is the need for jump or branch type microinstructions. 

After a certain amount of trial and error, the microinstruction 
word was divided into three parts. The first part consists of five 
bits which specify a basSic micro-opcode. The next three bits are used 
as an opcode modifier. The last part consists of eight bits which may 
specify the operand address or an eight bit operand itself. 

To increase the overall speed of the computer, the idea of per- 
mitting concurrent operations was considered. By studying the basic 
matrix operations, it was observed that quite frequently more than one 
of the address registers required incrementing between successive data 
transfer or arithmetic steps. A Similar requirement was found to exist 
for the counters. Therefore a variation of the microinstruction format 
was developed using each bit in the last two portions of the word to 


specify concurrent operations. To illustrate this idea, an example 
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using the "add One Direct" (AOD) microinstruction is shown. 
Example: 


AOD 0000 0 
u 


0 
(not used) K 


1 0011 
J AX YZ 
This microinstruction will add one to the contents 
of the J, X, and Z registers simultaneously. 
A companion microinstruction, "Subtract One Direct'' (SOD), was also 
established. 

A third concurrent operation which appeared useful allows the 
Simultaneous transfer of the contents of the A register to several 
different locations. This permits the referencing of this quantity 
at later points in the microprogram. The microinstruction is "Store A" 
(STA) and has a similar format to the AOD instruction. 

The remaining opcodes which were established are of a sStraight- 
forward nature and appear in Appendix III with detailed instructions 
for their uSe. 

Microprogram Structure 

This third and final part of the Computer Design concerns the 
actual structure and timing of the complete microprogram necessary to 
perform a main instruction. Looking at the main program, it appears 
that each main instruction is performed using the two-cycle system of 
RNI and EXEC described earlier. Actually this requires two separate 
microprograms. The first microprogram performs the RNI cycle. It trans- 
fers the first word! of the main instruction to the R register, further 
transfers the two address portions to specific locations, and then re- 
peats the whole operation for the second word of the main instruction. 


lrach main instruction usesS two consecutive words in the main 
memory. 
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Actually the same RNI microprogram is used for every instruction which 
means that a given part of every main instruction is always transferred 
to the same location during the RNI cycle. A detailed description of 
the RNI microprogram is shown in Appendix IV in a symbolic form using 
the notation found in ee 

The EXEC cycle is performed by a separate microprogram which of 
course depends on the operation required. The main opcode actually 
"calls" the appropriate microprograms as mentioned earlier. 

Looking now at the microprogram structure, it is found that each 
microprogram consists of a sequence of microinstructions and in general 
appears very similar to programs used in modern general-purpose compu- 
ters. Here again the basic two cycle system is used but in this case 
the cycles are performed by means of hardwired logic as in the conven- 
tional computer. The actual writing of the microprograms is exactly 
the same aS writing a conventional subroutine. The entering arguments 
which in this case are the operand addresses, are in known locations 
and the desired location of the result has also been specified. 

To illustrate some of the concepts and procedures that have been 
developed in this section, an example is presented in Appendix IV. 

The main opcode calls for a transfer of a block of data from the Main 
Memory to the Control Memory. The starting address of the block in M, 
the number of words in the block and the starting address of the block 
in N are also specified in the main instruction word. The micropro- 
gram RNI cycle is described and is understood to occur between every 


microinstruction. 
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4. Computer Test 


Now that a microprogrammed computer has been designed, it is 
necessary to apply a test to evaluate the results. Of course, one 
could determine the absolute times required for various operations, 
but this would not indicate the true worth of the concept. A more 
meaningful procedure would be to compare the microprogrammed computer 
against one of conventional design, intended for similar usage. 

The computer chosen for this comparison was the UNIVAC 1830 Avion- 
ics computer. This computer 1s a miniaturized, solid state, general- 
purpose computer, advertised for use in airborne and missile command 
and control systems. This of course, is the type of usage envisioned 
for the microprogrammed computer and is the chief reason for choosing 
the 1830 for comparison tests. The 1830 1s an outgrowth from the NTDS 
Unit Computer, AN/USQ-20, which was also built by UNIVAC and it uses 
the same instruction set. A complete technical deScription and list- 
ing of the instructions can be found in [3]. The basic characteristics 
and instructions have been excerpted from [3] and are presented in 
Appendix VI. 

It 1s common practice when comparing one computer versus another 
to prescribe certain "benchmark" programs. These "benchmark" programs 
are special tasks which occur frequently in computer usage and require 
more than the simplest programming to implement. They usually involve 
a certain amount of "bookkeeping" and arithmetic computation. In gen- 
eral, they are chosen to exercise as many of the hardware and software 
characteristics of a computer as possible. One of the most frequently 
used "benchmark" programs is matrix multiplication. This not only 


checks the efficiency of the multiplication routine, but also shows 
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how well a computer can do multiple addressing and loop indexing. 

Each "benchmark" program is programmed in the language of the 
computers to be compared, and then the running times are either calcu- 
lated or measured. The resulting times are usually a good basis of 
comparison. Other factors, such as ease of programming or flexibility 
of instructions, may also be important depending on the user's appli- 
cation. 

In this particular case, where the computers are intended as ele- 
ments of a control system, matrix arithmetic appeared to be a source of 
appropriate "benchmark" programs, since many of the control laws are 
formulated as recursive matrix equations. 

As a preliminary test, the five baSic matrix operations (addition, 
subtraction, transposition, multiplication, and inversion) were pro- 
grammed for both computers. These programs will handle matrices up to 
the dimensions of 7 x 7 and are listed in detail in Appendix V. 

The major "benchmark" program chosen for the test was the Kalman 
Filter. The filter is a set of three recursive matrix equations, 
which predict the state variables of a system based on noisy measure- 
ments of one or more of the variables. The particular feature of this 
filter is its ability to reduce the variance of the predictions due to 
a noisy environment. A detailed description of the filter can be found 
in Pe Appendix VII lists the equations, explains the notation used 
and the equivalent notation that appears in the programming, and gives 
the matrix dimensions assumed for the calculation of running time. 

The relative timing required for both the preliminary programs 
and the major benchmark" test is presented in the Evaluation and 


Conclusions section. 
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5. Evaluation and Conclusions 

The basis for evaluating the microprogrammed computer waS a compar- 
ison of the execution times required by the microprogrammed computer and 
the 1830 computer for each test program. For each microprogram, a for- 
mula was derived which gives the execution time in terms of minor cycles 
as a function of the operand matrix dimensions. Similar formulas were 
derived for each subroutine with the results in terms of main memory 
cycles. These formulas are shown in Table I. The variables, n and o, 
are subscripted in the formula for MATML since this is the only program 
that can have two operands of different dimensions. 

The first or preliminary test was to compare the execution times 
for each of the five basic matrix operations. All operand matrices 
were assumed to be 4 x 4 and using the formulas mentioned above, these 
times were calculated and are also shown in Table I. 

The final and major test was to compare the execution times for 
the "benchmark" program, the Kalman Filter. Using the matrix dimen- 
sions shown in Appendix VIII, the time to complete one iteration of 
the filter equations was computed for each computer and the results 
appear as the final entry in Table I. 

Up to this point no absolute time scale had been assigned to the 
microprogrammed computer. However, to make the comparison of the test 
results more meaningful, it was necessary to do so. Since the adver- 
tised memory cycle time for the 1830 computer was 4 microseconds, it 
seemed only fair to assign that Same value to the microprogrammed 
computer's main memory or major cycle. This then sets the minor cycle 
at .4 microsecond, based on the earlier assumption that the control 


memory was to be an order of magnitude faster than the main memory. 
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These timing assumptions were used to convert the test results in cycles 
to absolute values of time, which appear as the other two data columns 
in Table I. 

The first observation one might make is that the microprogrammed 
computer appears to be 3-5 times faster than the 1830 computer. This 
was the anticipated result but the magnitude of the difference should 
not be considered as a general figure. It can be noted by looking at 
the subroutines written for the 1830 computer, that a considerable nun- 
ber of instructions were used to transfer the arguments into the Sub- 
routine. This was necessary since the 1830 has no provisions for in- 
direct addressing. As a result the subroutine execution times are 
somewhat longer than might be expected on a machine having this feature. 
The percentage of the total execution time consumed by this task will 
decrease however, as the matrix dimensions increase, Since the argu- 
ment location is performed only once, whereas the computational in- 
Structions are repeated depending on the matrix dimensions. 

A second meaningful observation, though perhaps not apparent, 
iS a comparison of the number of cycles required for each program. 
These figures as shown are indicative of the number of instructions 
or steps required by each computer to perform the given matrix opera- 
tion. It will be noted that the microprogrammed computer requires 
approximately 3 times the steps or instructions to accomplish the same 
task as the 1830. This fact bears out the contention of this author 
that effective utilization of the microprogramming concept requires the 
use of a high-speed control memory. Quite obviously, the micropro- 
grammed computer would have been approximately 3 times slower than the 


1830 if the microprograms had been stored in the main memory, or an 


31 


auxiliary memory of the same speed. 

Several conclusions are apparent at this point, one of which is 
that the memory sizes should be increased. The control memory was not 
large enough as originally designed to contain the five matrix micro- 
programs, much less allow for temporary storage areas. Also the main 
memory size of 256 appears small aS was mentioned earlier. It is felt 
that making the main memory 4096 words and the control memory 1024 words 
would be adequate. Using these two figures, several other changes are 
also recommended. 

The main memory word should be increased to 20 bits. With 12 bits 
required to address 4096 words, 8 bits remain for opcodes. This 8 bit 
section could be further divided to allow 6 bits to specify 64 entry 
points in the control memory and 2 bits as the opcode modifier. The 
ability to transfer matrix elements into the control memory, when per- 
forming an operation requiring multiple accessing of a given element 
Still seems desirable, and thus the control memory word should be at 
least 20 bits in this case. 

Actually there is a trend today by those manufacturers utilizing 
microprogramming in some form to use a long control word. The advan- 
tage of this is to make possible a greater combination of concurrent 
operations specified by one microinstruction. Looking back at the 
microprograms, one can see that there are many sequential operations 
performed which are independent of each other, and thus could occur 
simultaneously if the appropriate instruction could be formulated. 

Tf a 20 bit control memory were used, 10 bits of course would be 
required to address any location in the control memory. The remaining 


10 bits could be subdivided to give 6 bits of opcode specification and 
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4 bits as an index or indices. This would allow 64 opcodes which should 
be adequate to designate all the necessary operations for general pro- 
gramming as well as those needed for I/O and special microsubroutines 
such as multiply and divide. 

The use of the dual-purpose address registers, X, Y, and Z, and 
the counters, I, J, and K, along with their initialization registers 
was very convenient and should be retained. The X, Y, and Z registers 
would have to be increased to 12 bits in this proposed version. 

Another question, regarding the advantages of microprogramming, 
concerns the amount of logic circuitry required. Generally speaking, 
the microprogramming concept is supposed to essentially substitute 
software for hardware. Of course, eliminating the operation counter 
would cause this reduction, but a certain amount of logic is necessary 
to implement the specified microinstructions, particularly those calling 
for various combinations of concurrent operations. 

Whether there would be a net reduction of logic in the micropro- 
grammed computer as designed cannot be answered definitely without first 
doing a complete circuit design. 

One final consideration is that of size and weight. The advertised 
figures for the 1830 are 1.1 cubic feet and 50 pounds. It is felt that 
the microprogrammed computer should be no larger than this, and with 
the present advances in integrated circuits, would probably be even 
smaller in physical size. 

In conclusion, it is felt that implementation of the micropro- 
gramming concept of computer design can provide a user with a special- 
purpose computer capable of being reprogrammed to fit a variety of 


applications at a competitive cost. 
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TABLE I 


Microprogram Execution Time Formulas 
MATAD 
112 + 43nm 
MATSB 
112 + 41nm 
MATPS 


128 + 40mn + 4m - 12n 


229 Pins [1 7n, + 4 + mp(80np + 19)] + 3np [14mg - 4 | + 4m, 


MATIN 


to 


120n° + 116n2 + 21n + 62 


Subroutine Execution Time Formulas 
MATAD 

78 + 11inm 
MATSB 

78 + 11nm 
MATPS 

ZO Maier al See So 
MATML 

75 + ng [18 + mp(3ing + 25)| 
MATIN 


39n3 + 91n? + 12n + 104 
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Computer Test Results 
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Operation Microprogrammed Computer UNIVAC 1830 Computer 
minor cycles microseconds major cycles microseconds 

MATAD 688 275 254 1016 
MATSB 656 262 254 1016 
MATPS 7 36 294 435 1740 
MATML 8981 5592 Zoo 10124 
MATIN 14646 5858 4104 16416 
Kalman 
PLlcer 40977 16391 13142 52568 
Program 
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APPENDIX I 


COMPUTER BLOCK DIAGRAM 


This Appendix contains a block diagram showing the physical 


organization of the computer and the data and address channels. 
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APPENDIX II 
MAIN INSTRUCTION FORMAT 
This Appendix contains a description of the main instruction for- 
mat and a listing of the main instruction opcodes. 
A main instruction actually consists of -two sixteen-bit computer 


words. The division of these two words into fields is shown below. 





Me tt "d tt 
————————>_>---—- 
First Word Ayes [7 ]2]s oft psa papwayis [6] 
Ha" a5 b 
—_——~ a 


i a 
ells lsl|s[shoha 2 papepsfie 


The 'ta'' and "i'' fields contain the operation code. These codes 


Second Word 





are described below. The ''b'' field contains the address of the first 
operand and the ''c'"' field contains the address of the second operand, 
if any. Finally the "d"' field contains the address at which the re- 

sult is to be stored. 

Although each main opcode appears as an eight bit quantity, it 
actually consists of an "a'' and an "i'' field. The five bit "a" field 
specifies the microprogram entry point and the "1'' field is used as a 
branching indicator within a microprogram. 


The following is a list of the main opcodes used in this project. 





Opcode Field 
iat oe 
MATAD 00101 001 
MATSB 00101 010 
MATPS 00110 001 
MATML 00111 010 
MATIN 01000 000 
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APPENDIX IIT 
MICROINSTRUCTION FORMAT 

This Appendix contains a description of the microinstruction for- 
mat and a listing of the microinstruction opcodes. 

Shown below are the two possible formats for a microinstruction 
word. The "ff", "j'', "y'" configuration is used where the instruction 
entails a single operation. The "f", "k'"' configuration is used with 
instructions calling for concurrent operations. 


hae oe a an 
oo OrOV——. 


1 2]3}4}5{6{7} [9 ]}rof1a Jz fr] 4]is [6 
ee ee k3ky ks ke k7 kg kg ky oki1 


The "k'' field is actually composed of eleven single bit fields 






Microinstruction Word 


labelled kj through kj}. 

The micro-opcodes used in this project appear on the following 
pages with, the details for their use. The actual bit patterns are 
shown in octal to clarify the presentation. 

LDA j y Ol Load A 
This instruction transfers the sixteen bit quantity found in the 


location specified by j and y to the A register. 


j = 1 andy =1: Load (M<x>) 
\ eee Load (M<Y?) 
y = 3: Load (M¢Z>) 
y = 5: Load (N<CX?) 
¥o=—20: Load (NCY>) 
varus: Load (NCZ>) 
y = ll: Load (X) 
y = 12: Load (Y) 
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y = V3: Load (Z) 


y = 15: Load (S) 
y = 16: Load (Q) 
) = 23 Load the contents of cell in 
N whose address is given in y 
j= <3? Load the quantity in y 


STN j y 02 Store N 
This instruction transfers the sixteen bit quantity in the A 
register to the cell in N addressed by y. The j field is not used 


in this instruction. 


SAL j y 14. Store A Lower 
This instruction transfers the contents of the lower half of the 
A register to the lower half of the cell in N specified by y. The j 


a -* 


field is not used. 


ADD j y 03 Add to A 
This instruction adds to the A register the quantity in the loca- 
tion specified by j and y. The codes for the j and y fields are the 


same as for the LDA instruction. 


SUB j y O4 Subtract from A 
This instruction subtracts from the A register the quantity in 
the location specified by j and y. The codes for the j and y fields 


are the same as for the LDA instruction. 


IPU= jy 12 Unconditional Jump 
This instruction causes a unconditional jump to the microinstruc- 


tion in the cell addressed by y. The j field is not used. 
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APPENDIX IIT 
MICROINSTRUCTION FORMAT 
This Appendix contains a description of the microinstruction for- 
mat and a listing of the microinstruction opcodes. 
Shown below are the two possible formats for a microinstruction 
word. The "ff", "j", "y'" configuration is used where the instruction 
entails a Single operation. The "f", "k'" configuration is used with 


instructions calling for concurrent operations. 


Fa i lal a 
1 [2]sJaTs| 6[o[s[s hola fra fasfaafas fie. 
ee k3ky ks ke k7 kg Kg ky okiy 


The "k'' field is actually composed of eleven single bit fields 





Microinstruction Word 


labelled kj through kj,- 

The micro-opcodes used in this project appear on the following 
pages with, the details for their use. The actual bit patterns are 
shown in octal to clarify the presentation. 

LDA j y 01 Load A 
This instruction transfers the sixteen bit quantity found in the 


location specified by j and y to the A register. 


7 = 1 and y = 1: Load (M<x>) 
y = 2: Load (M<Y>) 
y = 3: Load (M<Z>) 
y = 5: Load (N<X>) 
zee Load (N€Y>) 
y=7: Load (N<Z>) 
y=1l: Load (X) 
y = 12: Load (Y) 
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y = 13: Load (Z) 


y = 15: Load (S) 
y = 16: Load (Q) 
a = 22 Load the contents of cell in 
N whose address is given in y 
jcaeos Load the quantity in y 


STN j y 02 Store N 
This instruction transfers the sixteen bit quantity in the A 
register to the cell in N addressed by y. The j field is not used 


in this instruction. 


SAL j y 14 Store A Lower 
This instruction transfers the contents of the lower half of the 
A register to the lower half of the cell in N specified by y. The j 


field is not used. 


ADD j y 03 Add to A 
This instruction adds to the A register the quantity in the loca- 
tion specified by j and y. The codes for the j and y fields are the 


same as for the LDA instruction. 


SUB j y 04 Subtract from A 
This instruction subtracts from the A register the quantity in 
the location specified by j and y. The codes for the j and y fields 


are the same as for the LDA instruction. 


JPU J y 12 Unconditional Jump 
This instruction causes a unconditional jump to the microinstruc- 


tion in the cell addressed by y. The j field is not used. 
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ZIP jy ll zero Jump 
This instruction causes a jump to the microinstruction in the 
cell addressed by y if the counter specified by j is zero. Otherwise 


the next microinstruction is executed. 


j = 1: Sense K, jump if zero 
ae Sense J, jump if zero 
ia: Sense I, jump 1f zero 


NIP jy 13 Non-Zero Jump 
This instruction causes a jump to the instruction in the cell 
addressed by y if the counter specified by j is non-zero. Otherwise 


the next instruction is executed. 


j=l: Sense K, jump if non-zero 
j = 28 Sense J, jump if non-zero 
j = 4: Sense I, jump if non-zero 


AJP jy i A Jump 
This instruction causes a jump to the instruction addressed by y 
if the conditions of j are satisfied. Otherwise, the next instruction 


is executed. 


j = 1: Jump if A = 0 
j = 2: Jump if A # 0 


TP le yyy 10 Index Jump 
This instruction causes a jump to the instruction in the cell 
addressed by y if the contents of j field is identical to the i field 


of the main instruction contained in the R register. Otherwise the 
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next instruction is executed. 

The next group of micro-opcodes are those which specify concurrent 
operations. Each element of the k field specifies a particular opera- 
tion and by setting any combination of the k field elements to 1, any 


desired group of concurrent operations can be specified. 


STA k 10 Store A 

This instruction transfers the sixteen bit quantity in the A 
register to the locations indicated by the k field elements which are 
set to l. 


ee Store A in M<¢X> 
only one element may be set 


Ko = 1: Store A in M<Y> 
at a time 
k 3 = ls: Store A in M<Z> 


ieee ee Store A in NX> 
only one element may be set 


ks = 1: Store A in N<Y> 
at a time 
ke = 1: Store A in NCZ 
k7 = 1: Store A in X 
ko = ls: Store A in Y 
kg = l: Store A in Z 
k = ls: Store A in S 
k = 1: Store A in Q 


Note that only one element of kj, ky, and k3 may be set to 1 at a time. 


This is also true for the ky, k5, and k¢ group. 


AOD k 05 Add One Direct 
This instruction increments by one the contents of the register 


indicated by the elements of the k field which are set to l. 
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l 
Ko 
not used 
K 
Ky) 
ks = ls: Increment [I 
k 6 = 1; Increment J 


K7 = 1 Increment K 
Ke = 1: Increment A 
Kg = 1: Increment X 
Kigoe 2s) lnerement. \ 


k = l: Increment Z 


SOD k 06 Subtract One Direct 
This instruction decrements by one the contents of the registers 
indicated by the elements of the k field which are Set to 1. The k 


field element codes are the same as for the AOD instruction. 


RXF k 07 Right Transfer 

This instruction transfers the quantity in an initialization 
register to its companion counter or address register. Any combina- 
tion of the six possible transfers may be Specified by setting the 
appropriate k elements to l. 

Example 
= 1 and k = 1 


k 


6 re 


This means transfer the contents of ngs to Y and the contents 


of K, to K Simultaneously. 





i 
A 
Ky 
ky, = 0: 
ee ee 
ke = 1: 
aes 
ke = 0: 
kg = is 
Kio = i 
Ki =e 


not used 


(always zero) 
Transfer (X,) 
Transfer OE®, 
Transfer (2) 
(always zero) 
Transfer (T,) 
Transfer (J,) 


Transfer (K,) 


to 


to 


to 


to 


EO 


to 
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APPENDIX IV 
A MICROPROGRAMMING EXAMPLE 

This Appendix shows an example of the microprograms required to 
execute the RNI and EXEC cycles for one main instruction. The timing 
of the similar program cycles is treated in detail and the notation 
used to present these programs is explained. 

The time required to perform a given operation in a computer is 
generally dependent on the memory access time, therefore it is common 
to specify operation time in terms of the memory access or cycle time. 
In this computer there are two memories and it has been assumed that 
the Control Memory cycle time is one-tenth that of the Main Memory. 
To differentiate between the two different cycle times, the terms major 
cycle and minor cycle are used. One major cycle equals ten minor 
cycles. Thus to read a cell in M requires a time equal to ten minor 
cycles or, aS more commonly stated, ten minor cycles. Likewise to 
read a cell in N requires one minor cycle. 

To perform a microinstruction, it must be first transferred from 
N to the decoding register F and at this same time the Control Memory 
Address Register, H, is incremented by one, preparatory to reading 
the next instruction. This is the microprogram RNI cycle, which is 
implemented by hardwired logic. This cycle is understood to occur 


between each microinstruction and is shown symbolically below. 


(NCH?) —}F, (H) + 1—>H. 
This, be definition, requires one minor cycle. The time required to 
decode and execute a microinstruction is considered to require one 
minor cycle provided it doesn't involve a transfer to or from Main 


Memory. If this is required, then the microinstruction is considered 
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to require one major cycle. 

In the following example and all subsequent program listings the 
Operation Time given is in terms of minor cycles and includes the RNI 
cycle. 

For the purpose of the example, it is assumed that a block of 
data is to be transferred from a location in M to a location in N. 

The main instruction opcode for this operation will be TRNFR. The ''b" 
field will contain the Starting address of the block in M, the ''c" 
field will give the number of words in the block, and the ‘''d" field 
will specify the starting address in N. 

The first microprogram shown will perform the RNI cycle. This 
program is the same one intended for use in the computer and will be 
stored in the first three cells (000-002) of Control Memory. The pro- 
gram is listed in terms of symbolic statements since no formal mnemonics 
are associated with these microinstructions. All operations listed in 
the same cell occur simultaneously. 

AS an aid to understanding the symbolic notation, the following 


list shows the computer status at the end of this (RNI) microprogram. 


Register Contents 
C Address of next instruction 
F Last microinstruction of the RNI 
microprogram 
H Address of the entry point for the 


TRNFR microprogram 


R Marts... ls. andy Da £i1elds, GE tie main 
instruction 

X The starting address in M of the data 
block 

Y The number of words in the data block 

Z The starting address in N of the data 
block 
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The TRNFR microprogram performs the main program EXEC cycle. 
First the number of words in the data block is put in counter I, then 
a loop is used to transfer each word via the A register. Each time 
through the loop, the address registers are incremented, and the counter 
decremented and sensed for zero. At the end of the transfer, the pro- 
gram jumps to cell #000 to commence the next main program RNI cycle. 

All numbers shown in the program listing other than the times are 


in octal. 
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Oper Cell 
Time # 

11 000 

11 001 

2 002 

Oper Gell 
Time a 
2 15 

2 100 

A 101 

11 1G2 

11 103 

2 104 

2 105 

2 106 

2 107 


RNI MICROPROGRAM 


Contents 


(Mnemonic) 


Contents 


(Symbo lic) 


(M<C>)——> RR, (C) + 1-——C 


(A,.--Ag)——-JY, (Ag---Ayg)—2Z, 


1 
(C) + 1—>C, (M&C?)—IR 
(Ay---A5)—pH, (Ag---A 6) 9X, 


(Cc) +1— $C 


TRNFR MICROPROGRAM 


Contents 


JPU O 100 


LDA 1 012 


SIN GL 


LDA 1 O01 


STA 0 040 


AOD 0 006 


SOD 0 100 


NJP 4 302 


JPU 0 000 


(Mnemonic) 
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Contents 
(Symbolic) 


(100) —>H 

Ow) 

(A) — I 

(MX?) —— A 

(A) —NCZ? 

(X) + 1— SX, (Z) +152 
(I) - 1-31 

102——)H 


I # 0: 


0—H 


APPENDIX V 
MATRIX ARITHMETIC PROGRAMS 

This Appendix contains a listing of the programs written for the 
microprogrammed computer and the 1830 computer which allow both com- 
puters to perform the same basic matrix operations. Throughout the 
program listings, a matrix operation is denoted by the generalized 
algebraic equations or transformations shown below. A, B, and C de- 
note any matrix with the necessary conformability aSsumed, and their 
use aS operands is understood to mean the address of the first cell of 


the memory block containing the matrix elements. 


MAT AD A+ Be=C 
MATSB A =B = ¢ 
MATPS Aa—pal 
MATML AxBe=C 
MATIN A— a7! 


The notation used in presenting the microprograms for the micro- 
programmed computer is the same as in the example program in Appendix 
IV. It will be noted that there is only one microprogram for MATAD 
and MATSB. This was possible since the "bookkeeping" was the same for 
both operations, and only the arithmetic operation was different. The 
microprogram in general is a straight-forward implementation of the 
following algorithm. 

City 

The MATPS microprogram is also straight-forward and uses the 
algorithn, Bij = a where B = Al. 

The MATML operation was accomplished in a slightly different 


manner than usual. Since each element in a matrix iS referenced more 
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than once in the course of a matrix multiplication, it was decided to 
transfer both operand matrices into the Control Memory prior to begin- 
ning the multiplication, since the access time there is one-tenth that 
of the Main Memory. Further, it was found advantageous to transpose 
the B matrix while in the course of transferring it to N. Therefore 
the MATML microprogram uses a portion of the MATPS microprogram as a 
microsubroutire. 

The matrix inversion microprogram also transfers the operand to 
N, augmenting it in the process, and then uses the Gauss-Jordan elimi- 
nation method to accomplish the inversion. Finally, the elements of 
the inverted matrix are transferred back to M. 

To provide the necessary multiply and divide capabilities, two 
microsubroutines were written using the "shift and add" and "shift 
and subtract" methods respectively. The starting addresses of these 
microsubroutines are assumed to be #300 for Multiply and 4330 for 
Divide. 

It will be noted that all microprograms are shown as starting 
in cell #100 of the control Memory. Of course, each program would 
normally be located in its own section of Control Memory, but it was 
found at this point that the memory, as originally designed, was not 
large enough to contain all the microprograms. Rather than redesign 
the entire computer, the fact 1S recognized and the microprograms are 
presented in this general manner. 

The programs for use with the 1830 are presented in the format of 
subroutines. This was done to maintain the Similarity with the micro- 
programs, and generally these programs are available in this form for 


most general-purpose computers. Also included at the beginning of each 
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subroutine, are the instructions which must appear in the main program 
to "call" the subroutine. The general method of performing the var- 
ious matrix operations is the same as used in the microprograms with 
the exception of matrix multiplication. Here the usual algorithm 
stated earlier is implemented directly. The use of parentheses in the 
mnemonic form of insStructions indicates those operand values which 
must be inserted by the subroutine itself based on the entering argu- 
ments. 

Preceding each microprogram and subroutine is a brief flowchart 


to assist the interpretation and understanding of each program. 
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. ___ MATAD+ MATSB Nieeopragam — 





MATAD and MATSB Microprogram 


Oper Cell Contents Contents 
_Time_ _# _(Mnemonic) _ __(Symbolic) 
ial 100 LDA 1 001 (M<X?) ——A 
11 101 STA O 402 (A) —>M<Z>, (A)—OS 
2 102 AOD 0 007 (X) + 1— 3X, (Y) + 1-9 Y, (Z)+192Z 
BL 103 LDA 1 002 (M<Y>) ——>A 
sl 104 STA 0 401 (A)—>Q, (A) —IM@ 
2 105 LDA 3 060 110— A 
2 106 SAL 0 303 (Ay) ——? 303, 
56 (Multiply microsubrout ine) 
2 107 JPU O 300 300-—>>H 
2 110 Sito i (A) -—J I 
2 ae AOD 0 007 (X)+1—9X, (Y)+1—Y, (Z)+1—92Z 
11 52 LDA 1 001 (MCX?) —DA 
2 113 JPI 2 066 (i) = 2; 116— >H 
isl 114 ADD 1 002 (A) + (M<Y)) — >A 
2 nS: JPU O 067 117—9 4 
io 116 SUB 1 002 (A) - (M<Y?) —A 
vi ky STA 0 400 (A) —M<Z? 
2 120 SOD 0 040 (J)-1- 949 
2 1a NJF 2 111 3 #0: 11l1— 5H 
2 122 JPU 0 000 000— > H 
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a = i" 
oe 


Calevlate 
total 


term 


STrone 
in J 
© 
N= Oe 
Store N 6 Store 
ee 


&) 


VY 


S 


Ox) 


Return toled 
MATINL 
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Oper Cell 

Time # 
Z 100 
2 101 
2 102 
Tt 103 
11 104 
2 105 
2 106 
2 107 
ll 110 
2 ey 
i a2 
2 1s 
2 114 
2 115 
Z 116 
2 ey 
56 
Z 720 
2 ulzaa| 
2 22 
2 bps: 
2 124 
11 125 
Z 126 


MATPS Microprogram 


Contents 


LDA 1 013 


STN 


AOD 


LDA 


STA 


STN 


AOD 


RXF 


LDA 


STN 


STA 


AOD 


AOD 


LDA 


SAL 


3 


0 


Z 
fe) 
007 
001 


402 


004 


020 


001 


W 


401 


001 


005 


095 


303 


JPU O 300 


(Mnemonic) 


Contents 
(Symbo lic) 


CA Beer 

Ch) 7 oe 

(X) + 15X, (Y)4+1(9Y, (Z)+1-92 
(MCX?) — A 

(A)—3S, (A) MZ) 

(AY —JK, 

(X) + 1— x 

(Z)—7Z 

(M{X)}-—A 

(= 

(A)——}Q, (A)—JM(Z) 

C7) + leap, 

(X) + 1—— 5X, (Z) + 1— z 
120——jA 

A; —> 303, 

300-— H 


(Multiply microsubroutine) 


SIN O° J 


LDA 


STN 


RAL 


1 


0 


0 


SOD 0 


LDA 


JEL 


I 


2 


O11 
Xo 

001 
020 
001 


106 


57 


Oo ar a) 
(X) —— >A 

CO =a Xe 
(K,)—> K 
(K) ~ 1—>K 
(MKX) )——>A 


Gi)=2s) TS la—25H 


Oper Cell Contents Contents 


Time # (Mnemonic) (Symbolic) 
11 127 STA 0 400 (A)——> M¢{Z? 
2 130 JPU O 107 132—H 
11 131 STA 0 O40 (A) — NZ) 
2 132 SOD 0 O40 (J) - 1-3 J 
2 1533 WIE 22 125 J = 0: 150— 3H 
2 134 ZJP 1 117 K = 0: 142——>H 
2 135 AOD 0 001 (Z) + 1—2 
2 136 LDA 1 011 (X)——> A 
2 137 ADD 3 W (A) + (W)-— 9A 
2 140 STA 0 020 (A)—— > x 
2 141 JPU O 101 124— H 
2 142 LDA 2 X, (X,)— OA 
2 143 AOD 0 101 (A) + 1—— A 
2 144 STN 0 X, (A)—3X, 
2 145 RXF O 100 (XK 
2 146 AOD 0 001 (Z) + 1——9Z 
2 147 JPU O 100 123———H 
2 150 JPI 10 (i) = 1: O——7H 
2 sal JPU O Gig en 


58 





Cro pro 


Np XM p 


Mires o 


Nax Me 


1B] 


Thana 
Ape im N 


Sf 


59 





MATML Microprogram 


Oper Cell 
Time # 
tl 101 
11. 102 
2 103 
11 104 
Z 105 
2 106 
2 1G 
11 110 
ap rat 
Z IL 
2 Ls 
2 114 
2 115 
56 
2 116 
i rie 
2 120 
2 it 
2 122 
2 123 
56 
2 124 
2 125 
2 126 
Z 27 


Contents Contents 
(Mnemonic) (Symbolic) 

LDA 1 001 (MCX) )—> A 

STA 0 400 (A) —>M¢Z) 

STN O I, Qa 

LDA 1 002 (M6 = A 

STN O K, CO Fae JI 

STA 0 002 (A)— 5S 

AOD O 007 (X) + 1— 5X, (CY) +1— Y, (Z)+1 52 
LDA 1 002 (M{Y))——>A 

STA 0 401 (A) —>M@), (A)- 30 
STN 0 W (A)—> W 

LDA 3 145 116— A 

SAL 0 303 (A, )—9303,, 

JPU 0 300 300—> H 


(Mulitply microsubroutine) 


SEN Ore (A)—7J 

Oo 
LOA 21 OO er a 
STA O OO1 (A)—» 0 
LDA 3 153 124——A 
SAL 0 303 (A, ) 303, 
JPU O 300 300-5 


(Multiply microsubrout ine) 


STN O I = 
LDA 213 (Z)—O A 
STN O Z, 0 er 8 
LDA 3 TEMP TEMP —A 


61 


Oper 
Time 


Cell 


130 


spew 


es2 


Vee, 


134 


136 


137 


140 


141 


142 


143 


144 


145 


146 


147 


150 


jes 


12 


133 


154 


155 


156 


157 


160 


Contents 


(Mnemonic) 


STA 0 004 


STN 


AOD 


LDA 


SOD 


AOD 


NJP 


SAL 


JPU 


0 


0 


) 


0 


0 


ik 
006 
001 
100 
060 
162 
013 
Yo 


174 


MATPS 
151 


100 


( ) 


Contents 


(Symbolic) 
(A) zz 


(A) —T 

(X) + 9X, (Y) + 1-9 Y 
(MCX) )——>A 

(I) - 1 I 

(X) + 1—X, (2) + 1—42 
I # 0: 133——H 

Ci) mer a 

Oy. 

145 —A 

(A, )—? (MATPS 151) 


MATPS 
( 123 )— H 


(Transpose microsubroutine) 


LDA 


STN 


LDA 


STN 


RXF 


AOD 


AOD 


LDA 


STA 


LDA 


2 


0 


A 
Xo 


000 


164 


002 


001 


001 


006 


005 


002 


006 


62 


(T)— A 

CO? ke 

Omaar 

(A)—>T 

(Xo )—9X, (Y,) PY, (Z,) 92 


(I5)— 1 
C008 rome) 


Ce a 

(Z) + 1—z 

CX) 4K Ye Lae 
(N€X))— A 

Os 


NC) 


Oper 


Time 


Cell 


161 


162 


E62 


164 


165 


166 


L67 


L270 


17 


72 


eS 


174 


| 2s, 


176 


LFA 


200 


201 


202 


203 


204 


205 


Contents 


(Mnemonic ) 


STA O 001 


IDA 3 194 


SAL0 303 


JPU 0 300 


(Multiply 


ADD 


STN 


SOD 


NJP 


STA 


ZIP 


SUB 


STA 


SOD 


ZIE 


SUB 


STA 


JPU 


2 


0 


T 


al 


060 


184 


400 


207 


O11 


020 


182 


100 


000 


O12 


010 


Ee 


Contents 
(Symbolic) 


(A)— > Q 
69> A 


(A, ) —?303L 
300 —>H 


subroutine ) 


A) Crk 
(A)——>T 

(K) - 15K, (J) - + >J 
K #0: 155——H 
> a 

(A) —9 M¢Z2) 

J = 0: 200—5H 
OO 

Ce aa 
(A)—5X 

15 3——>H 

(I) - 1—JI 

I =O: O—-+H 
(Qs 

(A) IS 
O20 er 4 
152—H 












| 
| 
| 
= aa 
J,°7 an | 
W =, |= TEMP | 
a ! 
: 
| | 
BAD tS: 
Xo= 2 P 
J=To5 Y=Y% } 
| 
r : | 





LL 


Yar J X *Xo 
K=Ko9 S290 





@ 


Y= TEMP 
VFAN 


E 
Zz. rae 
Q NI 


EX 
X=X+Z -7 


K=2 ? 


ie sVe+Z 


=< 
i 
ox 





MATIN MICROPROGRAM 





Oper Geil Contents Contents 
Time # (Mnemonic ) (Symbolic) 
v1 100 LDA 1 001 (M<X?)——-A 
2 101 SIN. O. 1 (A) =e 
2 102 SINO K, (A) —>K, 
11 103 STA 0 400 (A) —— MZ) 
2 104 ADD 2 I, (A) + (1,)— A 
2 105 STINO J, (A) — 5 
2 106 SIN O T =—7 
2 107 AOD 0 003 (X) + l— xX, (Z) + b—9Z 
11 110 LDA 1 001 CM 7k 
1 111 STA 0 400 (A) —M<Z) 
2 ee LDA A013 (Z)— A 
2 113 STN O Z, 29 ae 1 
2 ee LDA 3 001 1——~A 
2 115 STA 0 004 (A) ——>Z 
2 116 LDA 3 TEMP TEMP — > A 
2 ipl? STINO wW (A)——W 
2 1X0) SNOW. C7 oe 
2 121 SOD 0 010 (A) - 1—JA 
2 1? STA 0 010 ()—— yy 
2 123 RXF 0 001 (K,)—>K 
2 124 RXF O 006 Go eG 
2 125 AOD 0 006 (X) + +—->xX, (Y) + 1-7 Y 
iy 126 LDA 1 001 (MAX? )——-A 
D 127 STA © 100 (A) —> N¢YD 


6/ 


Oper Cell Contents Contents 





Time a (Mnemonic) (Symbolic) 
2 130 SOD 0 140 Gi) = lees (I) = lee 
2 131 NIP 4 134 I #0: 125——H 
2 132 LDA 2 K (K)— >A 
2 leg SUBE eeu (A) - (J)— OA 
2 134 AJP 1 137 A=0 : 137—H 
2 135 LDA 3 000 O— >A 
2 136 JPU O 140 140—H 
2 137 LDA 3 001 1——A 
2 140 AOD 0 002 (Y) +1 Y 
2 141 STA 0 100 A— NY) 

2 142 SOD 0 O40 (J) - 1 3I 
2 143 NIP 2 132 JI #0: 132—H 
2 144 SOD 0 020 (K) - 1—— 3K 
2 145 NIP 1 124 K #0: 124-—H 
2 146 [DAN 2 Ya Cae 
2 147 STN 0 Xy C= 
2 150 RXF 0 046 (YY, (Ig) J, (Ip JI 
2 151 LDA 1 006 (N<Y))—F A 
2 152 STAs 0 802 (A)— >S 
2 153 LDA 1 006 (NCY)) A 
2 154 LDA 3 143 143——A 
2 155 SAL 0 333 (Ay, ) 9 3331, 
2 156 JPU O 330 330—H 
68 (Divide microsubrout ine) 
2 157 LDA 1 016 (0) —A 


68 





Oper 
Time 





162 


163 


164 


oS 


200 


201 


202 


203 


204 


205 


206 


207 


Contents 


(Mnemonic) 


STA 0 100 


SOD 


AOD 


NJP 


RXF 


LDA 


SUB 


SOD 


AP 


LDA 


STA 


LDA 


STA 


LDA 


STA 


LDA 


SAL 


0 


0 


3 


0 


040 


002 


153 


143 


020 


174 


O11 


aL 


020 


005 


002 


006 


O01 


203 


302 


JPUOr 300 


(Multiply 


STA 0 001 


LDA 1 005 


SUB 1 016 


STA 0 200 


SOD 0 040 


69 


Contents 


(Symbolic) 
aN 
(Ce er) 

(Y) + 1— 9Y 
JI 40: 153— 5H 


OS ae SF ee SEY ae. K, 
(X,) —aX 


(K) ——>A 
(A) - (I)-—— SA 
(K) - 1——>K 


A #0: 174—H 
OO mer 
CeCe 
== 
(N<X)) —A 
(A) ——S 
(N<Y>)——>A 
(A)— 0 
203—PA 

(Ay )—— 3031, 
300-— H 


subroutine ) 


(A) —Q 

N ((X>)—— A 
(A) - (Q)—— >A 
(A)——9N<xX) 
(609 Mee reg 2) 


Z 210 
2 Zi) 
2 242 
2 23 
2 214 
Z 215 
Z 216 
2 217 
fp 220 
2 221 
2 222 
2 225 
2 224 
2 225 
2 226 
2 22) 
Z 230 
Z 231 
Z 232 
2 235 
2 2 34 
2 235 
2 236 
2 237 
Z 240 


Contents 


(Mnemonic ) 


AOD 0 006 


NJP 


LDA 


SOD 


AJP 


JPU 


SOD 


ZIP 


ADD 


STN 


RXF 


LDA 


AOD 


STN 


SOD 


STN 


JPU 


2 


1 


176 


013 


010 


O11 


020 


O01 


225 


C12 


165 


100 


241 


OA 


013 


002 


(Real 


70 


Contents 


(Symbolic) 
OO) +) le, CY) I 


J #0: 176— JH 
2) 

(A) - 1——7#A 
Coes Os 
(A) —>x 

1— 3A 

(A) - (K)}— 3A 
A=0 :;: 223——H 
(Yo) -—HY, (J9)—I 
165 —H 

(I) - 1—JI 

I = 0 : 241— 5H 
Oa 
(7 
Cae 

0 aaa ak 

Oy ay 31 

Ay, 4A A, 2) eZ 
9 cay 2,05 

OF ras ts 

(A) - 1——>A 

(A) WJ, 

PP oan) 
151—H 








Oper Cell 

Time # 
2 241 
2 242 
Z 243 
2 244 
2 245 
2 246 
2 247 
2 250 
2 255 
2 252 
2 253 
2 254 
2 255 
2 256 
2 257 
2 260 
2 261 





Contents 


LDA 


STN 


STA 


LDA 


STA 


AOD 


SOD 


NJP 


SOD 


LJP 


2 


0 


(Mnemonic ) 


020 


005 


020 


005 


O40 


250 


020 


000 


002 


O11 


246 


71 


Contents 
(Symbolic) 


Cl) aA 

(A)—>J, 

CLG? oy CK) a a ee 
(Z) + 1—>Z 

(Ww) —>A 

(A) + (I,)—— pA 

Aas 

(N<X))—A 

(A) ——M (Z) 

(X) + 1—aX, (Z) + 1—35Z 
(J) -1—3J 

d 4-0: 250-—pu 

(K) - 1—K 

K=0 : 0O— >H 

UE) er 5) 

(X) ——>A 

246 ——>H 


MULTIPLY Microsubrout ine 


Oper Cell 
Time os Contents (Symbolic) 
2 300 C= ©, 167-76 
Oia 0s 20 loa (G) 9 ce 
ne bt CAT: 20, .) KAS = --045) 
Q16 = 1: 302—>H, (G) - 1— 6G, 
G=0: 303— OH 
2 302 (A, ---Ayg) + CS) (An. - - -Q))5 
(Q1---Qy5)-—9(Qy-+ O16), 301—H 
2 30 3 ( ae tl 


Average Operation Time: 44 minor cycles 
Maximum Operation Time: 68 minor cycles 

This Multiply microsubroutine is composed of special-purpose 
instructions which, like the RNI microprogram, have no mnemonic 
names assigned. This method of implementing a multiplication 
is necessary to make the multiply time competitive with conven- 
tional machines. The G refers to a special 3-bit counter used 
with the multiply and divide routines exclusively. 

The first instruction performs the necessary initialization. 
The second instruction accomplishes all the steps which must be 
performed each time, i.e., sensing the right most bit in the Q 
register for 0 or 1, and decrementing the counter. If 6 0, 
the contents of AQ are shifted right one bit and the same instruc- 
tion is read and executed again. The reading and executing take 


two minor cycles. 


72 





eae 1, the next instruction is executed also requiring 2 
minor cycles. This instruction adds the multiplicand in S to the A 
register, shifting the result one bit to the right, and at the same 
time shifts Q one bit right. 

Then a jump back to #301 occurs and the cycle is repeated until 
the multiplication is complete. At that time instruction in cell #301 
senses G = 0, and causes a jump to #303. Cell #303 contains a standard 
JPU instruction whose operand is set by the calling microprogram and 
provides the return to that microprogram. Note that the total execu- 
tion time depends on the number being multiplied. An average and a 
maximum time are shown. A value of 56 minor cycles was used in cal- 


culating program running times. 


ne 


DIVIDE Microsubroutine 





Oper Cell 
Time a Contents (Symbolic) 
2 330 0-7 0216 ==2 G 
2 Sal (Ay---Q)) SCS) Anite A ie (G eas | 
A>0: 10),, 3314 
GFO0: 
APO: (A)+(S)\—S, 00, ,, 
2 SoZ 331—Y H 
G=0:; 333-——H 
2 S35 ( —- H 


Operation Time: 68 minor cycles 

This Divide microsubroutine is also composed of special-purpose 
instructions to permit competitive timing. 

Again the first instruction performs the initialization. The 
second instruction performs the equivalent of left shifting AQ one 
bit and then subtracting S from A all in one operation. The third 
instruction senses the contents of A and takes the appropriate action 
to implement the "shift and subtract" procedure. This two instruc- 
tion cycle is repeated until the division is complete. The G 
counter is sensed for 0, and when found, a Jump to the final in- 
struction occurs. This last instruction is a standard JPU which 


causes a return to the calling microprogram. 
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MATAD + MATSB Svbr opting 





MATAD and MATSB Subroutines 


Oper Cell Contents Contents 
Time FF (Mnemonic) (Symbo lic) 
3 600 RTJ 000 100 *600"-—9100, , 101—P 
= 601 1 000 2000 ho = 1000-ne= 2000 
: 602 2 - -"1 3000 +1 for Add, -1 for Sub, C=3000 
= 100 0 002 (600) 
3 101 RPL.Y+1 010 100 (100,) + 1— 100, 
2 102 ENT.A 010 °#100 (100; )——>a 
2 103 STR.A 0:10 114 (AJ —1 14, 
2 104 STR.A 010 #116 (A) — 116, 
S 105 RPL.Y+1 010 100 GC0 Da a7 100", 
2 106 ENT.A 010 #100 (100; )—>A 
2 107 STR.A 0 10 (Sari 
2 110 STR.A 0 0g (A) 123, 
3 ial RPL.Y#1 010 100 (100,) + 1——>100, 
2 112 ENT.A 010 #4100 (100; )——A 
2 13 STR.A 01°0°. 133 (A) — 133, 
2 114 ENT.A 0 2-07 — 601) (60114) —A 
2 1s STR.A Oily Ge 140 (A) ——>140, 
2 116 ENT.A 010 (601) (6017 )——3 A 
2 7 STR.A O:a70 144 (A)—141, 
2 120 STR.A 010 142 (A) — 142, 
2 121 ENT.A 020 (602) (6021) —PA 
2 n22 STR.A 010 137 (A) —A1 371, 


76 


Oper 
Time 


Cell 
# 





123 


124 


125 


127 


130 


Tesi 


gZ 


133 


134 


135 


L3G 


Voy 


140 


141 


142 


143 


144 


145 


146 


Contents 

(Mnemonic) 
ENT.A 010 
STR.A OrEG 
ENT.B" 200 
ENT .Q Gres 
STR.Q 0 32 
BSK.B" 2 00 
ENT.A 032 
STR.A O32 
MUL O32 
ADD.Q 000 
STR.Q 010 
BSK.B™® 20 0 
ENT .Q 000 
ENT.A 0.0 6 
ADD.A Deo? 
SUB.A Oi23-2 
STR.A O22 
BSkKvn =20020 
JMP 12080 
JMP ior 0 


(602) 
143 
0000 
(1000) 
(3000) 
0 
(1000) 
(3000) 
(1000) 
2 
144 


0000 


(+1 or -1) 


(1000) 
(2000) 
(2000) 


(3000) 


(m x n+2) 


140 


(60 3) 


7] 


Contents 
(Symbolic) 


(6027) a 

GQ) aarlts 

0000—>B7 

(1000)——Q 

(Q) —3000 

(B2) + 1—>B2 

(1001)——>A 

(a2 3001 

(A) x (1001)——>AQ 

(Q) + 2-0 

(Q) —— 144, 

(B2) + 1—>B? 

+1 or -l1—40 

(1002)—>A, Skip NI if Q Neg 
(A)+(2002)—A, Skip NI if Q 
Pos 

(A) - (2002)— A 
(A) — 3000 


Skip NI if B2 =mxn + 2, 
or (B2) + 1—>B2 


Go to 140 


Return to Main Program 









(a0c0+B*) = 
(1000+8*) 


(B7= B+ (k52) 
BY B*+7 


D) 


(155) = (5) +7 


B= (sas 
(153) =(153)- Z 


(Q-——* (153) =O iu 


"e 


VY 


MATPS Subroutine. 








Oper Cell Contents Contents 
Time £ (Mnemonic ) (Symbolic) 
3 600 RTI 000m) 100 "600" S100, & “10 e 
z 601 1 000 2000 #£=A = 1000, A! = 2000 
. 100 0 000 (600) 
3 101 RPL.Y+t1 010 = 100 (100,) + 1— 100, 
2 102 ENT .A 0. 10 LOO (100, )——>A 
2 103 STR.A 010 £110 (A)——>110 
2 104 STR.A 01-0 es OS rl UC 
3 105 RPL.Y+1 010 100 (100,) + 1——5100, 
2 106 ENT.A 010 100 (100, )——>A 
2 107 STR.A 010 151 (A)—151, 
g) 110 ENT .A 020 (601) (601),)——>aA 
2 La STR.A 010 # 121 QO 2, 
2 112 STR.A OC: 0 wees (= =7l35, 
2 113 ENT.A 010 (601) (601,)——>A 
2 114 STR.A 010 £122 (A) —122,, 
2 1S STR.A 010 # £136 (A) —1 36, 
2 116 ENT .Q 010 =| -1—~>Q 
3 117 ENT.B" 100 0000 o— >s! 
3 120 ENT.B? 200 0001 I—> BZ 
ist Pass 2nd Pass 
2 121 ENT .A 0 31 1000 (1000)——pA (1001)——A 
2 ia) STR.A 363-2 2000 (A) —}2001 (A) 2000 
2 123 STR.A DESO 1152 "skap" (A) ——} 152 
2 124 STR.A 000 0000 (A) ——0 "skip" 
3 125 BSK.B” 100 0005 (B4y+41— 8 (ply+1—>B! 
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Oper 
Time 


5 


Cell 
# 





126 


27 
P30 
131 
132 
133 
134 
135 
136 
137 
140 
141 
142 
143 
144 
145 
146 
147 
150 
es | 
152 
13 


154 
sys) 


Contents 


(Mnemonic) 
BJP.B 2 OO 
ENEGE = =29000 
ENT.A 000 
STR.A 030 
ENT.A Oegn0 
STR.A Ois300 
STR «0 030 
ENT.A OF 391 
STR.A 032 
STR Bos “a7 <0 
ADD .A OmeE 
ENT.B" 170 
BSK.B" 200 
RPle vol) Ona. 0 
JP 500 
RPL.Y+l 0 30 
ENT.B" 130 
RPL.Y-1 0 30 
JP 500 
JP 100 


121 


Z 
0002 
155 
152 
153 
154 
(1000) 
(2000) 
0000 
2 
0000 
0000 
154 
135 
iS 
2 
I Ses: 
134 


(602) 


Contents 
(Symbolic) 

lst Pass 
(B2)-1—>B2 
jump to 121 


en 

A 

(A) 9155 
Gs2)==s7 

Co =S153 

(Q) ——7 154 

(1000 + B!)——»a 
(A) —> 2000 + B2 
(B!)—-yA 

CA) eel 2) —<— 
(Se 

(Gy Ss B- 
(154) - 1+—— 154 
Go to 135 if (A) 
(155) + 1— 3155 
Gis ==> Be 
(153) - 1— 153 


Go to 134 1£ (A) 


Go to 602 


Temporary Storage and Counters 
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2nd Pass 
go to NI 





hon zero 


non- Zero 


So a ee — LLL SS LS A A Se a A pS 


_.  LMATML Subre otin a 





(a0v)=Me@ 
(Q00)= RX 


eB? = (20 6) 
B= (208), (203g 
2 


= 

Co 
=D 
@_ 
i> if 
x i, 
Q> 

&> 

Cc 

oS 
&, 
Xo 








(3000+ 8*)= 
(2000+8°)+Q | 


B= B'+z 
B*=B°+Ms 
£3) =(203)-2 


~, 








(QnE\V ALAS jth 
(179) = (177)-2 


a > eee et 
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Oper Cell 
2 600 
- 601 
- 602 
- 100 
] 101 
2 102 
2 103 
S. 104 
2 105 
2 106 
2 107 
2 ae 
2 111 
2 2 
2 1S 
2 114 
2 115 
2 116 
2 iD, 
2 120 
2 2A 


RTJ 


0 


RE Ga 


STR.A 


SER eet 


REE. Y 


SRA 


ENT .A 


SER. 


SERA 


ENT .A 


oO TR.A 


SiR A 


SlRoe 


ENT .A 


STR.A 


SLRs 


SLR ek 


STR.A 


MATMUL Subroutines 


Contents 
(Mnemonic) 
07070 
000 
000 
000 
+1 010 
O° 140 
dae o 
+1 010 
010 
O22 40 
Gale} 
010 
010 
010 
010 
oa hee 8, 
9 ip fee. 8 
010 
010 
010 
Oe 8) 


100 
2000 


3000 


(600) 
100 
106 
111 
100 
115 

(601) 
124 
155 

(601) 
130 
ile is) 
156 
602 
125 
134 
151 


160 


84 


Contents 
(Symbolic) 


"6000"—}100, '101"—>P 


A 


1000, B = 2000 


C 3000 


(100,) + 1-100, 
(A) — 106 
(A)——7111 

(100,) + 1—7100, 
(A) ——y115 
(601y)——>A 

(A) ——J1 24 
(A) — 155 

(601, )——jA 
(Aj=-—7, F350 

(A) — 5133 
(A)——156 

(602; )—>A 
(A) — 125 

(A)——7 134 

(A)— 151 
(A= 60 





Oper 
Time 


3 


3 


Cell 


ee 


122 
igi 
124 
125 
126 
p27 
130 
131 
eo 
133 
134 
135 
136 
137 
140 
141 
142 
143 
144 
145 
146 
147 
150 
151 


eZ 


Contents 


(Mnemonic) 
ENT.B™ 250 
ENT.B" 20 
ENT .A 0 3 
STR.A Ons 
STR.A 0 3 
BSK.B™ 30 
ENT.A 0 3 
STR.A OG 
BSK.B" 20 
ENT.A 0 3 
STR.A 0 3 
STR.A 0 3 
ENT .A 00 
STR.A 0 3 
RPL.Y+1 01 
STR.A 01 
ENT.A 00 
STR.A 0 3 
ENT. A 0 3 
STR.A Gag 
ENT.B" m3 
BSK.B" shail 
ENT .A 01 
STR.A 0 3 
ENT .A 0 3 


O 0000 
0 0000 
O (1000) 
QO ¢ 3000) 
0 Le? 
O 0005 
2 (2000) 
0 201 
QO 0005 
24.2000) 
3 ¢€ 3000) 
0 202 
O 0002 
0 205 
QO 5000 
0 176 
QO 0002 
0 206 
O 202 

O 204 
QO 206 
O 0000 
QO 0000 
3 (€ 3000) 
0 201 
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Contents 
(Symbo lic) 


038° 
een 
(1000)-——_A 
(A) ——>3000 
(A)-—177 
(B>) + 1—»B? 
(2000)——->A 
(A) -— > 201 
Gye 
(2001) ——A 
(A) 3001 
(A) —>202 
2— 3A 
(A) —— 205 
(100,) + 1—100, 
(A) -—— 176 
2—— A 

(A) ——>206 
(202)—_54A 
(A)-——9204 
(206) B? 


(Bo) lee 


G=-- A 


(A)——3000 + B? 
Cie 


Oper Cell Contents Contents 


Time _# (Mnemo nic) (Symbolic) 
2 153 STR.A Ouse 203 (A) —> 203 
3 154 ENT .BD 130 205 £(205)— >B! 
2 155 ENT.Q 0 31 (1000) (1000 + B!)y—~40 
10 156 MUL 0 32 (2000) (2000 + B2) x (Q)——aQ 
- 157 RSH. AQ On) 0) 24, Right Shift AQ "z" bits 
3 160 RPL.Y+0 0 33 (3000) (Q)+(3000 + B3)——>3000 + B® 
3 161 BSK.B® 110 0000 Bly + 1—yB! 
3 162 STR.B™ 270 0000 (B2)—>A 
2 163 ADD.A 0-370 02 (A) + (202)——>A 
3 164 ENT.B” 270 0000 (A)—>p2 
3 165 RPL.Y-1 030° 9203 (203) - 1— 4203 
2 166 JP sy 0/08) By s Jump to 155 if (A) non zero 
3 167 RPL.Y#1 030 #£206 (206) + 1—206 
3 170 RPL.Y-1 0 8.0 W048 (204) - 1—3204 
2 Lal JP 500) 146 Jump to 146 if (A) non zero 
2 P72 ENT .A O30 201 (201) — A 
3 173 RPL.A*Y 030 205 (205) + (A)—— 205 
3 174 RPL.Y-1 03-04 “177 (177) - 1—>177 
2 LS JP 5 0 0 142 Jump to 142 if (A) non zero 
2 176 JP 100 (603) Jump to 603 
lye 
201 
202 Temporary storage 
203 and Counters 
204 
205 
206 
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MATIN Subrovtine 


_ — See 
aa 





_—— 


(270) =(40%0 +@7) 


(4000+B* ) = 


(4000+ @*) 220) 


Bt = 4 + 7 
far!) =@.2/)-2Z 


ao —* ~eeeeeeee = 8 ee ee ee ee 2 ere ee + re 0A ES LS 





(270) =(4000+8) 


¥oootR') = 
ihe _- | 
(2926) (4o00 +7) | 


B'=R'+2Z 
B>= B°+Z 
a7) = (87/7) =2 


\ 





2 
fA 


(R ) 88 





mmm em a ES RG A SS SS a 


(274 2. +7 
(069) = (847) -2 
26b) = (2665-2 


© <n 
ie 


Themahin A 2 


sy 2000 


MATIN Subroutine 


Contents 
(Symbolic) 


Bogie 100hen OI e 


1000, A7! = 


Oper Cell 
Time _# 
| 600 
= 601 
- 100 
3 101 
2 102 
2 103 
3 104 
2 105 
2 106 
2 107 
2 110 
2 Pie 
2 112 
2 13 
2 114 
3 115 
é' 116 
2 dial 
2 120 
3 21 
2 122 
Z T23 
2 124 


Contents 
(Mnemonic) 

RTJ 000 100 
1 000 2000 
0 000 #4600 

RPL.Y +1030 °#100 

STRoA 030 #106 

STR.A 62320" 

RPL.Y +1030 #100 

STR.A O<3 08265 

ENT.A OR 2 OrmarGeOle 

STR.A Oa Om 117 

STR.A O70 LS 

ENT .A 010 (601) 

STR.A 010 120 

STR.A Onle0: 1122 

STR.A 010 £256 

ENT .B™ 100 0002 

ENT.B™ 200 0000 

ENT.A 0 30 (1000) 

STR.A 0 3 2 (2000) 

BSK.B” 300 0005 

STR.A 0 3 2 (2000) 

STR.A 030 £266 

STR.A OS. Orme? 72 
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(100) + 15100 
(A) —— 106 
(A= 
(100) + 1—- 100 
(A) — 265 
(601,)-—>A 
(A) 117 

(A) —145 
(601;) ——3A 

(A) ——120 

CA) 7122 

0 ers aN 
2—>Bl 

0— B2 

(1000) —>A 
(A) —> 2000 
Ga ge ie 
(A) 7 2001 
(A) —> 266 


(A) 9272 


Oper Cell Contents Contents 


Time gue (Mnemonic) (Symbolic) 
2 125 ADD.A 030 266 (A) + (266)—— >A 
2 126 STR.A 030 267 (A) —> 267 
2 127 STR.A 030 277 (A) ——>277 
3 130 RPIy-1 0 3°60 272 (272) — 19272 
2 131 ENT.A 000 0000 O——A 
2 132 STR.A 030 2744 D277 
2 1633 ENT.Q 0 30 266 (266)—0Q 

10 134 MUL 0350 266 (Q) x (266)——>A0 
10 35 MUL 000 £0002 (Q) x (2)——>AQ 
2 136 ADD .Q 000 £0002 (0) + 2——0 

2 137 STR.Q Om 340 276 (OQ) 9.2.6 

3 140 ENT.B" 400 0002 — ne 

2 TEES ENT.A 030 266 (266)——>A 

2 142 STR.A 5356 DA G0 eee 0 

2 143 ENT.A 0-320 266 (266) —— A 

2 144 STR.A 030 270 (A) —— 270 

2 145  ENT.A 031 (1000) (1000 + Bly—~ya 
z 146  STR.A 034 4000 (A)—4000 + B4 
3 147. +=BSK.B2 100° 0000 (Bly + 1p! 

3 150 BSK.B2 400 0000 (Ge See 

3 ea REEs volo 0 270 (270) - 1—— 3270 
2 152 JP 5 0.0 145 A=0: 145— Pp 
2 13 ENT.A 0 30 266 (266)-—DA 

2 154 STR.A Ons 0) 270 (A) 5270 

2 155 ENT.A 030 270 (270)——yA 


ek 


Oper 
Time 


2 


2 


12 


Cell 


Sanaa 


156 


137 


160 


Tol 


162 


163 


164 


165 


166 


167 


170 


plea 


172 


bby: 


174 


173 


176 


177 


200 


201 


202 


20:3 


204 


205 


206 


Contents 

(Mnemonic) 
SUB.A 530 
ENT .Q 100 
ENT .Q 000 
STR.Q 034 
BSK.B" 400 
RPL.Y-1 0 30 
ape 5 00 
RPL.Y-1 0 30 
oe > 0210 
ENT .Q O70 
ADD .Q 000 
MUL O30 
RSH.AQ O 0 0 
ADD .Q 000 
ENT.B’ 10 0 
SER. 03-0 
ENT .A Oras o0 
STR.A 0 53°10 
ENT .A Ops 0 
SIR.A O30) 
ENT.A Osi 
STR.A F300 
ENT.A 02050 
ENT.Q 031 
DIV Oo an0 


oa gl 


0001 


0000 


4000 


0000 


270 


139 


27 


143 


207, 


0001 


274 


0000 


0000 


275 


267 


Ze 


272 


273 


4000 


270 


0000 


4000 


270 


92 


Contents 
(Symbolic) 


(A)-(271)—JA, Skip NI if A # 0 
1—5Q, Skip NI 

0-—9 

(Q)— 4000 

(eee 

(270) =a 270 


A #0 


155—>P 
1—9271 
143——5P 
(277) —Y0 

0) 1 


(271) 


A #0 


(Q) x (274)——AQ 
Right Shift AQ 
(Q) + 2— 
o—>B} 
(Nas 72 75 
(267)—A 
(A)-—— 271 

(272) —A 

(A) ——927 3 

(4000 + BL)—+A 
(A) —}270 
Cx 

(4000 + BLy—40 
AQ + (270) —3Q 


Oper 
Time 


2 


= 


Cell 


ae 


207 


210. 


211 


22 


Ns 


214 


2 


216 


Ze, 


220 


221 


22 


223 


22h, 


225 


226 


227 


2 30 


Zou 


232 


2353 


234 


239 


2 36 


2 oy, 


240 


Contents 

(Mnemonic) 
STR.Q 0 3 
BSK.B™ 10 
RPL.Y-1 0 3 
JP 5 0 
STR.B" iy 
ADD.A 0 3 
ENT .B® Ay: 
ENT .A 0 3 
STR.A Q 3 
ENT .A 0 3 
STR.A 0 3 
ENT.B" 3 3 
ENT.Q Ores 
MUL 0 3 
RSH.AQ 00 
RPL.Y-Q 0 3 
BSK.B"™ 10 
BSK .B® an 
RPL. Yo sv0N3 
JP 5 
STR=B 14 
SUB.A 53 
ENT.B" 10 
RPL.Y-1 0 3 
JP 50 
RPL.Yt1 O 3 


1 


4000 


0000 


271 


204 


0000 


274 


0000 


4000 


270 


267 


27 


209 


4000 


210 


4000 


0000 


0000 


Zr 


Za3 


0000 


276 


0002 


Pa es, 


213 


274 


a3 


Contents 
(Symbolic) 


(0) 4000 + BI 
(Bly) + 1—>pl 
(271) - 1-271 
A #0: 204-—P 
(Bly) ——a 

(A) + (274)-—A 
(A)——>B! 

(4000 + BLYN—> A 
(70 
(267)——>A 

029) oer cee 

(275) —B? 

(4000 + B3)—-4Q 
(Q) x (270)—+AQ 
Right Shift AQ 
(4000 + BL) - (9)—>(4000+B!) 
(BL) + 1—~ 5s! 
(B3) + 1—>B? 
(27) =. 2 271. 
A #0: 223—3P 
(B1)—yA 
G@=O76)— A skip NE oes 
Be 

(273) - 1— 9273 
A#0: 213—)P 


(274) + 1— 4274 


Oper Cell 


Time # 
S 241 
3 GN) 
2 243 
3 244 
3 245 
wy. 24.6 
3 247 
3 250 
2 25 1. 
3 252 
2 253 
2 254 
2 255 
2 256 
3 257 
3 260 
3 Zor 
2 262 
3 263 
2 264 
2 265 
266 
267 
270 
271 
272 
273 
278 
275 
276 
2/17 


Contents 
(Mnemonic ) 


RPL.Y-1 
RPL.Y-1 
JP 
BSK.B" 
RPL.Y*e1 
STR.A 
ENT.B™ 
STR.B"™ 
ADD .A 
ENT. B"™ 
ENT.A 
STR.A 
ENT.A 
STR.A 
Boks” 
BSK.B® 
RPL.Y-1 
JP 
RPL.Y-1 
JP 


JP 


0 


0 


326 


a0) 


267 
266 
167 
0000 
2H 2 
266 
0002 
0000 
266 
0000 
266 
273 
4000 
(2000 ) 
0000 
0000 
23 
255 
Z72 
250 


(602) 


Contents 


(Symbolic) 


67) a) 267 
(266) - 1— 266 
A #0 :@7—> P 
(B2) + 1——p2 
(272) + 1-—— 3272 
(A) ——}266 

2— B+ 
(Bt)—> A 

(A) + (266)— 3A 
(A) ——>B4 

(266) —>A 

(A) —727 3 

(4000 + B4)—aA 
(A)——>2000 + B2 
(Bo) + 1—>B4 
(B2) + 1—>B? 
(273) - 1— 273 


A #0 


255——> P 
(272) -— 1-272 
A #0 : 250—>P 
602——5P 


Temporary storage and Counters 
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APPENDIX VI 


TECHNICAL CHARACTERISTICS OF UNIVAC 1830 COMPUTER 


This Appendix contains a brief listing of the technical character- 
istics taken from [1], a description of the instruction word format, and 
a listing of the instructions used with the 1830 computer. 

Technical Characteristics 

Memory 
4096 words, 30 bits 
4 microsecond cycle time 

Control 
Single Address 
62 instructions 
Seven index register 
Seven branch designators 
Seven operand interpretation designators 


Instruction Execution Time 


Add: 8 microseconds 
Subtract: 8 microseconds 
Multiply: 32-48 microseconds 
Divide: 48 microseconds 
Jump: 8 microseconds 


Construction 


Semiconductor Microelectronic Integrated Circuits 
Instruction Word Format 


asf aa| ar [26] 25 en ]oa 2] ai [ao] fis 7[vopis] 


SS ee 
J k b 
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BERD PoPEP EEE PELL 


Function Code Designator 

Branch Condition DesSignator 
Operand Interpretation Designator 
Index Designator 


Operand Designator 
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APPENDIX VII 


KALMAN FILTER EQUATIONS 


This Appendix contains a listing of the Kalman Filter equations, 
the definitions of the terms, the equivalent terms used in the test 
program, and the dimensions used to compute the program running time. 
Equat ions 


* 
X(k/k) = X(k/e - 1) + CCR) [z(k) - Zek/e ~ 1] 
G(k) = | P(K/x -1)HT(HP(kK/k - 1)HT + R)~ "| 


P(k/k - 1) 


b [Pc Sh 2 pee Cee ene (ke Sy ee 2316" + 0 


1 


A * 
Z(k/k - 1) = HX(k/k - 1) = HGX(k - 1/k - 1) 


Definition of Terms 


Variable Definition Remarks 
$ nxn state transition constant for a given time- 
matrix invariant plant 
G nxl optimum filter gain 
matrix 
H plant observability matrix order of matrix depends 
upon number of plant 
observables 
R measurement noise covar- value selected by operator 


iance matrix 


Q random excitation dis- value selected by operator 
tribution covariance 
matrix 
E error covariance matrix initial value selected by 
operator 
Z nxl vector of noisy 


plant observable states 


99 


Variable Definition Remarks 
“A 
Xx nxl vector of the pre- 
dicted value of X(k), 
based on the previous 
observation X(k - 1) 


IN> 


nxl vector of the pre- 
dicted value of Z(k), 

based on the previous 

observation Z(k - 1) 


nxl vector of the best 
estimate of X(k), based 
on the current observa- 
tion of Z(k) 


[>< + 


Equivalence and Dimensions of Terms 


Programming 
Variable Equivalent Dimensions 

D PHI 4x 4 

G G 4x1 

H H 1x4 

R R Lee ee 

Q Q 4x 4 

1 le 4x4 

(hi Z se 

x XH 4x1 
A 

Z ZH. bE ccna 

* 

Xx XS 4x 1 
¢? PHIT 4x 4 
ue cei 4 xl 
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APPENDIX VIII 


KALMAN FILTER PROGRAMS 


This Appendix contains a listing of the Kalman Filter programs 
for both computers. These programs are a straight-forward implementa- 
tion of the equations shown in Appendix VII and in fact appear very 
Similar. The operation times shown are for the complete subroutine 
or microprogram "called", and are based on the matrix dimensions listed 
in the previous Appendix. 

Here again, Subroutine names and variables are used to denote the 
Starting addresses in memory of the appropriate programs or matrices. 
The terms TA, TB, and TC refer to sections of memory used as temporary 
storage. 

All constant matrices and the necessary transposes are considered 
to be resident in memory and the input matrix, Z, is also considered 


to be present when needed since I/0 programming was not considered. 
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Oper 


Time 


2095 


9005 


680 


9005 


9005 


712 


Zoe 


788 


179 


343 


2021 


for the Microprogrammed Computer 


Cell 
# 





00 


Ol 
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